Testing for dynamical dark energy models with redshift-space distortions 
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The red-shift space distortions in the galaxy power spectrum can be used to measure the growth 
rate of matter density perturbations Sm- For dynamical dark energy models in General Relativity 
we provide a convenient analytic formula of f{z)a8{z) written as a function of the redshift z, where 
/ = din (5m /din a (a is the cosmological scale factor) and as is the rms amplitude of over-density at 
the scale 8h~^ Mpc. Our formula can be applied to the models of imperfect fluids, quintessence, and 
k-essence, provided that the dark energy equation of state w does not vary significantly and that 
the sound speed is not much smaller than 1. We also place observational constraints on dark energy 
models of constant w and tracking quintessence from the recent data of red-shift space distortions. 



I. INTRODUCTION 



After the first discovery of cosmic acceleration from the distance measurements of the supernovae type la (SN la) 
[H, the existence of dark energy has been also supported from other observational data such as Cosmic Microwave 
Background (CMB) @, H and Baryon Acoustic Oscillations (BAO) [!]■ From the theoretical point of view, such a 
late-time cosmic acceleration is problematic because of the huge difference between the observed dark energy scale and 
the expected value of the vacuum energy appearing in particle physics Along with the cosmological constant A, 
many alternative acceleration mechanisms have been proposed, including modifications of the matter/energy content 
and large-scale modifications of gravity (see Refs. @ for reviews). 

The dark energy equation of state w is constrained by measuring the expansion rate of the Universe from the 
observations of SN la, CMB, and BAO Although it is possible to rule out some accelerating scenarios from 

the analysis of the cosmic expansion history alone, we require further precise observational data to clearly distinguish 
between models with subtly-varying w. So far the A-Cold-Dark-Matter (ACDM) model has been consistent with the 
data, but there are many other dynamical models such as quintessence |Tl|, k-essence [13: [3 1 and f{R) gravity 
which are compatible with current observations. 

The large-scale structure of the galaxy distribution provides additional important information to discriminate be- 
tween different dark energy models [15| . The galaxy clustering occurs due to the gravitational instability of primordial 
matter density perturbations. The growth rate of matter perturbations can be measured from the redshift-space dis- 
tortion (RSD) of clustering pattern of galaxies. This distortion is caused by the peculiar velocity of inward collapse 
motion of the large-scale structure, which is directly related to the growth rate of the matter density contrast Sm ila |. 
Hence the RSD measurement is very useful to constrain the cosmic growth history. 

Recent galaxy redshift surveys have provided bounds on the growth rate f{z) or f{z)as{z) in terms of the redshift 
z = 1/a — 1, where / = dln^m/^lna and erg is the rms amplitude of Sm at the comoving scale 8 Mpc {h is the 
normalized Hubble constant Hq = 100 /ikmsec"^ Mpc~^) |17H29| . Although the observational error bars of /erg are 
not yet small, the data are consistent with the prediction of the ACDM model [l^, . Recently the RSD data were 
used to place constraints on modified gravity models such as f{R) gravity and (extended) Galileons |30l. [Slj . Since 
the growth rate of matter perturbations in these models is different from that in the ACDM |33 - t34| . the allowed 
parameter space is quite tightly constrained even from current observations. 

For the models based on General Relativity (GR) without a direct coupling between dark energy and non-relativistic 
matter, the gravitational coupling appearing in the matter perturbation equation is equivalent to the Newton's 
gravitational constant, as long as the dark energy perturbation is negligible relative to the matter perturbation. 
Nonetheless the evolution of perturbations depends on the background cosmology, so that the dynamical dark energy 
models with w different from —1 can be distinguished from the ACDM. In particular, the future RSD observations 
may reach the level of discriminating between different dark energy models constructed in the framework of GR. 

In this paper we derive an analytic formula of f{z)as{z) valid for dynamical dark energy models including imperfect 
fluids, quintessence, and k-essence. Provided that the sound speed Cs is not much smaller than 1 and that the 
variation of w is not significant, our formula can reproduce the full numerical solutions with high accuracy. The 
derivation of f{z)as{z) is based on the expansion of w in terms of the dark energy density parameter fi^^, i.e.. 
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w = wq + ^n(f^2;)"- Since f(z)crs{z) is expressed in terms of the present values of and Vlx as well as w„ 

(n = 0, 1, 2, • • •), our formula is convenient to test for dynamical dark energy models with the observational data of 
the cosmic growth rate. For the models with constant w there are 3 par ameters wo, crs{z ~ 0), and f2i^(z ~ 0) in 
the analytic expression of ] (z)a^{z). In tracking quintessence models [35| , where the dark energy equation of state 
is nearly constant during the matter era (w = wjoj), we show that our formula of ]{z)g^(£) also contains only 3 
parameters: W(o), <y%(z = 0), and ^x{z = 0). Using the recent RSD data, we carry out the likelihood analysis by 
varying these 3 parameters to find observational bounds on wq and ui^o) . 

This paper is organized as follows. In Sec. |lT] we review cosmological perturbation theory in general dark energy 
models including imperfect fluids, quintessence, and k-essence. In Sec. Illll dvnamical dark energy models are classified 
depending on how the equation of state w is expanded in terms of ri^^. In Sec. IIVI we derive an analytic formula for 
/(z)t78(z) and in SeclVl^e confirm the validity of this formula in concrete examples of dark energy models. In Sec. lVIl 
we perform the likelihood analysis to test for constant w and tracking quintessence models with the recent RSD data. 
Sec. IVIIl is devoted to our main conclusions. 



II. COSMOLOGICAL PERTURBATIONS AND REDSHIFT-SPACE DISTORTIONS 

As the dark energy component we consider a fluid characterized by the equation of state w = Pxl Px, where 
is the pressure and px is the energy density. We also take into account non-relativistic matter (cold dark matter 
and baryons) with the energy density and treat it as a perfect fluid with the equation of state Wm = 0. We 
deal with such a two-fluid system in the framework of GR under the assurnption that dark energy is uncoupled to 
non-relativistic matter. Our analysis covers quintessence [llj and k-essence [l2l . [l3| models, in which the Lagrangian 
P of dark energy depends on a scalar field ^ and a kinetic term X = — (V0)^/2. In these models we have that P^ ^ P 
and px = 2XP^x - P, where = dP/dX. 

In the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) background with the scale factor a{t), where t is the 
cosmic time, dark energy and non-relativistic matter obey, respectively, the following continuity equations 

p'x + 3{l + w)px = 0, (!) 

Pm + 3Pm = , (2) 

where a prime represents a derivative with respect to A'^ = In a. We introduce the density parameters fl^ = 
SttGpx/ i'iH^) and ilm = 87rGp,„/(3if^), where G is the gravitational constant and H = a/a is the Hubble pa- 
rameter (a dot represents derivative with respect to t). From the Einstein equations it follows that 

nx + n,n^i, (3) 
^ = -lii + wnx). (4) 

The dark energy density parameter satisfles the differential equation 

n'^ ^ ^SwQxil - ^x) ■ (5) 

Let us consider scalar metric perturbations about the flat FLRW background. We neglect the contribution of tensor 
and vector perturbations. In the absence of the anisotropic stress the perturbed line element in the longitudinal gauge 
is given by [36j 

ds'^ = -(1 - 2$) dt^ + a2(t)(l + 2$)da;2 , (6) 

where <f> is the gravitational potential. We decompose the energy densities pi (where i — x, m) and the pressure Px 
into the background and inhomogeneous parts, as pi = pi{t) + 6pi{t,x) and Px = Px{t) + 6Px{t,x). We also define 
the following quantities 

. _ Spt „ _ V^w, 



Pi aH 



{i = x,m), (7) 



where Vx and Vm are the rotational-free velocity potentials of dark energy and non-relativistic matter, respectively. 
In Fourier space dark energy perturbations obey the following equations of motion (37l - l40l | 

d'^+3icl-w)Sx = -il + w){W + 9x) , (8) 
H' w' \ „ f kV f cl 
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where = 5Px/5px and fc is a comoving wave number. The perturbed equations of non-relativistic matter (perfect 
fluids) are 



5L 



3$' - , 



k 



From the Einstein equations we obtain 



k 



= - ( ^mSm + ^xSx I , 



$' + $ = -|f^) [n,ne,n + {1 + w)ax0x] 



where dm and Sx are the rest frame gauge-invariant density perturbations defined by 



aH 



= 4 + 3(^1 (l + w) 



(10) 
(11) 

(12) 
(13) 

(14) 



For imperfect fluids such as quintessence and k-essence there exist non-adiabatic entropy perturbations gene rated 
from dissipative processes. We introduce a gauge-invariant entropy perturbation Ss^ of dark energy, as (stI l38l. IZlj 



Px (5Px 5p, 



Px \ Px 



where Ca is the adiabatic sound speed defined by 



2 Px 

c„ = — = w - 



Px 



3(1 + If) ■ 



(15) 



(16) 



In the rest frame where the entropy perturbation is given by 5sx = (c^ — c^)5x, the sound speed squared = is 
gauge- invariant. Using the definition of 5x in Eq. ()14p . the pressure perturbation of dark energy can be expressed as 



5Px = cl Spx + 3 



V k 



{1+W)ici-Ci)pxi 



whereas the sound speed squared = SPx/Spx in the random frame is related to via 



In terms of cl, Eqs. © and ^ can be written as 



6'^ + 3(c2 - w)5x = -3(1 + w)$' - (1 + w) 



2 + ^-3cl]ex 



k 



2 / 2 
C 



1 + 9(X) 



(17) 

(18) 

(19) 
(20) 



l + w 

In k-essence characterized by the Lagrangian density P{(f),X), the perturbed quantities can be expressed as 

Spx - {Px + 2XPxx)SX^{P^~2XPx^)6(j), (21) 
SPx = Px(5X + P^J0, (22) 

&x = 



(23) 



where d4> is the field perturbation and 5X — (j)S4> + (fP^. From Eq. (|17l) the rest frame sound speed can be obtained 
by setting 50 = in Eqs. (PT|) and ((^ . i.e., = P^x/{P.x + '^XPxx) [Hi- In quintessence characterized by the 
Lagrangian P = X — V{4i), the sound speed squared reduces to = 1. 
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surveys 


0.067 


0.423± 0.055 


6dFGRS (2012) [26.1 


0.17 


0.51± 0.06 


2dFGRS (2004) [18] 


0.22 


0.42± 0.07 


WiggleZ (2011) [23.1 


0.25 


0.3512± 0.0583 


SDSS LRG (2011) [24.1 


0.37 


0.4602± 0.0378 


SDSS LRG (2011) [241 


0.41 


0.45± 0.04 


WiggleZ (2011) [231 


0.57 


0.415± 0.034 


BOSS CMASS (2012) [25] 


0.6 


0.43± 0.04 


WiggleZ (2011) [231 


0.78 


0.38± 0.04 


WiggleZ (2011) [23] 



TABLE L Data of /erg versus the redshift z measured from the RSD. 



From Eqs. (dU]) and (HH it follows that 

$" + (^2 + ) $' . (24) 



R ) I aH ) 



For the perturbations deep inside the Hubble radius (k ^ aH) relevant to large-scale structures, the r.h.s. of Eq. ([M)) 
can be neglected relative to the l.h.s. of it, in addition to the fact that (5m — ^m- If the contribution of dark energy 
perturbations can be neglected relative to that of matter perturbations in Eq. i.e. jr^m^ml ^ Eq. ([M)) 

reads 

Ka + ^ (1 - 3^"^-) ^'ra ^ ^ , (25) 

where we made use of Eq. 

During the deep matter era in which fim ~ 1, there is a growing- mode solution 5„i = (5^ oc a to Eq. (|25|) . In this 
regime Eq. ([T2|) tells us that $ = constant and hence Q^a — — 15' from Eq. (fTO|) . For the dark energy density contrast 
(5x , it is natural to choose the adiabatic initial condition [13, Eo] 

J, = (1 + w)b„, . (26) 

The initial condition of is known by substituting Eq. ((26)) and 8'^ ~ (1 + if + ti;')5m into Eq. (fT9|) . We will discuss 
the accuracy of the approximate equation (|25p in Sec. IIIII 

The growth rate of matter perturbations can be measured from the RSD in clustering pattern of galaxies because 
of radial peculiar velocities. The perturbation bg of galaxies is related to the matter perturbation a s bg = bbm, 
where 6 is a bias factor. The galaxy power spectrum Vg{k) in the redshift space can be expressed as jig. |43| - |45| 

- - '^^^^'Pge{k) + ^i^Vge{k) , (27) 

where jj, = k ■ r/{kr) is the cosine of the angle of the k vector to the line of sight (vector r), Vgg{k) and Vee{k) 
are the real space power spectra of galaxies and 0, respectively, and Vge{k) is the cross power spectrum of galaxy-0 
fluctuations in real space. 

In Eq. pO| the variation of the gravitational potential is neglected relative to the growth rate of bm, so that 

-fb,n , f=-^- (28) 

Under the continuity equation ([28)) . Vgg, Vgo, and Vee in Eq. (|27)) depend on (bbm)^, {bb,n){fbm), and (fbrn)^, 
respectively. We normalize the amplitude of bm at the scale 8h~^ Mpc, for which we write erg- Assuming that the 
growth of perturbations is scale- independent, the constraints on bbm and fbm translate into those on bag and fcr^. 
The advantage of using fas is that the growth rate is directly known without the bias factor b. In Table U we show 
the current data of /erg as a function of z from the RSD measurements. 



III. DYNAMICAL DARK ENERGY MODELS 



In this section, we discuss a number of dynamical dark energy models in which the field equations presented in the 
previous section can be applied. 
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FIG. 1: Evolution of fas versus the redshift z for = 1, fc = 0.07 h Mpc"\ asiz = 0) = 0.811, and = 0) = 0.73. 

The cases (a) and (b) correspond to w — —0.8 and w = —0.5, respectively. The solid lines show the numerically integrated 
solutions, whereas the bald dashed lines correspond to the solutions derived by using the approximate equation (|25|l . 



A. Imperfect fluids 



For imperfect fluids the rest frame sound speed Cs is generally different from the adiabatic sound speed Ca- For 
constant w one has = w. If is constant as weh, the evolution of dark energy perturbations is known by solving 
Eqs. and (PH)) together with the background equations ©-([S]). This is the approach taken in Ref. Note also 
that for Cs of the order of 1 the contribution of dark energy perturbations to $ in Eq. is negligibly small relative 
to matter perturbations [H, • 



In Fig. [T] we plot the evolution of /erg for 



1 and k = 0.07 h Mpc with two different values of w. The 



approximate equation (j25p reproduces the full numerical solution within the 0.1 % accuracy. This means that, for 
Cg = 1, the contribution of dark energy perturbations to the gravitational potential is suppressed relative to that of 
matter perturbations. For = and w = —0.8 we find that the difference of fos{z = 0) between the numerical result 
and the approximated solution is about 4% for the mode k = 0.07 ft. Mpc~^. For w larger than —0.8 the difference 
gets even larger, but such values of w are not allowed observationally 0, 0]- In Ref. [s^] it was found that galaxy 
redshift and tomographic redshift surveys can constrain the sound speed only if is sufficiently small, of the order 
of cl < 10-4. 



B. Quintessence 



Quintessence [Tl[ is characterized by the Lagrangian P(0, X) = X — V{(j)), where T^((/)) is the field potential. In this 
case the sound speed Cs is equivalent to 1, so that the contribution of dark energy perturbations to the gravitational 
potential is negligible. There are several classes of potentials which give rise to different evolution of w. 

The first class is the model with constant w, which can be realized by the following potential [l^ 



167rGf7^ 



■ sinh 



GttG 



4>Q + 0l) 



where a = (1 + 0o and fimo are the today's values of and fim respectively, and 



i + w 1 1 + yi - n„M 

6ttG \w\ \/^mO 

This case is identified as an imperfect fluid with = 1 and w — constant discussed in Sec. IIII Al 



(29) 



(30) 
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The second class consists of freezing models [48|, in which the evolution of the field gradually slows down at late 
times. A typical example is the inverse power-law potential [4^ 

= M^+P<j)-P , (31) 

where p is a positive constant. For this potential there exists a so-called tracker solution [s^ with a nearly constant 
field equation of state w = W(^q^ = —2/ {p+2) during the matter era, which is followed by the decrease of w. Considering 
a homogeneous perturbation Sw around i/;(o) , the field equation of state is expressed as [HO] 



w - u;(o) + ^ _ ^ ^ 2n(n + ( j ' ^^^^ 

Expansion of w around D.^ reads 



(l~w^o))u;(o) ^ (l-u;fo))u;fo)(8w(o)-l) 
2(1 - u.2ppi«3^)(4u.(0) - l)(18u-(o) + 1) 



nt + oi^t), (33) 



(1 - 2w^o) + 4.wf^))il - Sw^o) + 12u;2p))(l " Mo) + 24.wf^^) " 

which varies with the growth of fl^ . 

The third class consists of thawing models j48l | , in which the field is nearly frozen by a Hubble friction during the 
early cosmological epoch and it starts to evolve once the field mass m drops below H. In this case w is initially 
close to —1 and then w starts to grow at late times. The representative potential of this class is that of the pseudo 
Nambu-Goldstone boson [5l|, i.e., 

F(0) = [1 + cos(0//)] , (34) 

where A and / are constants. Assuming that the evolution of the scale factor can be approximated as that of the 
ACDM model, the field equation of state is estimated as [52| (see also Ref. [SS] ) 



w 



-l + (l + iZ;o)a'(^-i) 



{K - F{a)){F{a) + 1)^ + [K + F(a))(F(a) - 1) 



K 



_{K - n-J'\n-J" + i)A- + [K + n-J'\n-J" - 

where tDo and fJ^o are the today's values of w and respectively, and 



(35) 



(36) 



The constant K is related to the mass of the field at the initial displacement, (p ^ (pi. Expansion around ^Ix = ^ gives 
w = -l + wiVlx + \{K^ + l)winl + -^(3A"* + ilK^ + 3)wiO^ + 0{nl) , (37) 

O 1(0 

where 

The growth of VLx leads to the deviation from w — 

C. K-essence 

The equation of state for K-essence with the Lagrangian density P{(p,X) is given by w = P/{2XP^x — P)- This 
shows that cosmic acceleration with w k, ~\ can be realized either for (a) A sa or (b) Px ~ 0. 

In the case (a) Chiba et al. [l^l showed that, for the factorized function P{(j),X) = F(0)F(A), the field equation 
of state is given by the same form as Eq. ([55]) with the replacement K = -^/l — V^^{(l)i) / {6ttGF_x iO)V'^ {(pi)) , where 



7 



F{X) is expanded around X = as F{X) = Fifi) + Fx(0)X. In this case the sound speed squared is also close to 1, 
so that the situation is similar to that in thawing quintessence models. 

In the case (b) the evolution of w depends on the functional form of P{(f>,X), so it is difficult to derive a general 
expression of w [H^l- One of the examples which belongs to this class is the dilatonic ghost condensate model [H^ 
characterized by the Lagrangian P(0,X) = -X + e'^^'f'X'^/M'^, where k = VSttG, A and M arc constants (see also 
Ref. [15). In this model the fixed points during the radiation and matter eras correspond to Px — and w ~ —1, i.e., 
y = Xe'^^'^ /M^ = 1/2. On the other hand the (no-ghost) accelerated fixed point corresponds to = 1/2 + A^/(A)/16 
with the equation of state w = -[l- AV(A)/8]/[1 + i\^f{\)/S] > -1, where /(A) = 1 + y'l + 16/(3A2). Hence the 
evolution of w is similar to that in thawing quintessence models. 

The sound speed squared in the dilatonic ghost condensate model is given by cl = {2y — l)/(6y— 1), so that ~ 
during radiation and matter eras. The late-time cosmic acceleration occurs for 1/2 < y < 2/3 and hence < < 1/9 
at this fixed point. The fact that is close to during most of the cosmological epoch is a different signature relative 
to quintessence. However, since w is very close to —1 during radiation and matter eras, the adiabatic initial condition 
()26p shows that the dark energy perturbation Sx is initially suppressed relative to the matter perturbation 5m- As 
long as the today's value of w is not significantly away from —1 the contribution of the dark energy perturbation to 
$ is suppressed relative to the matter perturbation, so that the approximate equation (|25p can be trustable even in 
such cases. 



IV. ANALYTIC SOLUTIONS OF fag 

We derive analytic solutions of fas by solving the approximate equation (|25|) . We expand the dark energy equation 
of state in terms of the density parameter as 



Wo 



^w„(f^.)". (39) 



Since fix grows as large as 0.7 today, we expect that it may be necessary to pick up the terms higher than the first 
few terms in Eq. (j39|) . 

In terms of the function / = S'^/Sm, Eq. can be written as [57| 

Swnxil - ^-)^ = /' + ^(1 - SwQx)f - ^(1 - ^x) , (40) 
where we employed Eq. ([5]). Introducing the growth index 7 as / = {ilmV = (1 ~ ^x)'' [I^, Eq. (PU)) reads 

3wnx{i - nx) in(i -nx)-^ = ^- ^w{i - 27)17, + (1 - ax^ - ^(1 - a.)'"^ • (4i) 

We derive the solution of Eq. (pi]) by expanding 7 in terms of fix, i.e., 7 = 70 + J2n=i Ini^x)^- Using the expansion 
of w in Eq. p9p as well, we obtain the following approximate solution 

^ 3(1 -wq) 3 {1 - wo){2 ~ 3wo) + 2wi{5 ~ 6wo) 
^ 5-6wo 2 (5 - 6wo)^(5- 12wo) 

+ [{wo - l)(3wo - 2)(324u;g - 420wo + 97) + 12(5 - 6wo)^(5w2 - 12w2U'o + 12w?) 

+ 6(5 - 6u;o)(72u;^ - 90u;o + 23)wi]/[4(5 - 6wQf{5 - 12u;o)(5 - 18u;o)] + 0{nl) . (42) 

The 1-st order solution is identical to the one found in Ref. by setting wi = 0. For wq = —1, wi = 0, and W2 ~ 
it follows that 7 ~ 0.545 -I- 7.29 x 10~^i^x + 4.04 x lO^'^fi^. In this case the second and third terms are indeed much 
smaller than the first one, so that 7 is nearly constant. For the models with wq = —1, Wi = 0.3, and t«2 = (in which 
case the value of w today is around —0.8) one has 7 ~ 0.545 + 1.21 x 10~^ilx + 6.55 x 10~^f2^. Even in this case the 
variation of 7 induced by the second and third terms is small (see also Refs. [g^l for related works). 

From the definition of / the matter perturbation obeys the differential equation (InSm)' = {^ — ^x)'^- Using Eq. ([5]), 
we obtain 

'^15- ^xV~^ 
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In the following we derive the solution of this equation under the approximation that 7 is constant. We expand the 
term (1 — around fix = 0, as 



(1 - O^p-i = 1 + ^ a„(l),)". 



(-1)" 



(7-l)(7-2)...(7-n). 



(44) 



n=l 



In order to evaluate the r.h.s. of Eq. we expand l/w in the form 

00 



1 _ 1 

W Wo 



n=l 



(45) 



where the coefficients /3„'s can be expressed by Wn's, say /3i — —wi/wq. Then Eq. (|43)) can be written as 



l + ^c„(f2.)^ 



n=l 



i=0 



where 

n 

Cn — ^ ^ ^n—ii^i : 

with ofQ = /3o = 1- The first three coefficients c,; are 

Cl 
C2 



-(7-1) , 

Wo 

i(7-l)(7-2) + (7-l)-"^ 



Wo 



(46) 

(47) 

(48) 
(49) 



C3 



1. -.w „w „^ 1, -.w , , W2W0 - W3wg - 2wiW2Wo + , , 

-77(7 - 1)(7 - 2)(7 - 3) - -(7 - I)(7 - 2)— + (7 - 1) -2 -3 . (50) 



Integrating Eq. (|46l) . it follows that 



Wo 



(5,„o exp < - — 
3wo 



ln^ + ^^((l).o)"-(l^.)") 



17, 



(51) 



where (5^0 is the today's value of (Sm- Normalizing (J^o in terms of o^(z = 0), we obtain the following expression for 



/(T8 = (1 - ^xV <J%{z = 0) exp \ 

3wo 



(52) 



In terms of the redshift z the energy densities of non-relativistic matter and the field are given, respectively, by 
Pm ~ Pmo(l + and px = pxoexp[Jj' 3(1 + w)/(l + z)dz\. The 0-th order solution to the field energy density is 

obtained by substituting w = wq into the expression of px, i.e., pi^"* — Pxo{^ + z)^^^'^'"'°\ This gives the 0-th order 
solution to Qx, as fii"'' = r2a;o(l + z)'^^° /[I — flxo + ^^2:0(1 + 2)^'"'']- If we expand w up to first order with respect to 
fix, we can use the iterative solution w ~ wq + wiVl^x^ This process leads to the following integrated solution of fix ■ 

(53) 



^(1) ^ r!,o(l + ^)'"""' [1 - »xo + ^.0(1 + zf^o]^^/^o 

l~nxO + nxo{l + Z)3"'« [1 - flxO + nxoil + 2)3-o]-i/-o 



In the presence of the terms higher than second order, we can simply carry out the similar iterative processes. 

In tracking quintessence models the coefficients u>„'s (n > 1) are expressed in terms of Wo — W(o), see Eq. (|33p . 
Hence there are 3 free parameters cr^^z = 0), flxo, and wo in the expression of fag. In thawing models of quintessence 
one has u>o = —1 and Wn's (n > 2) can be expressed in terms of wi and K, see Eqs. (|37)) . Then fag in Eq. (|52]) is 
written as a function of z with 4 free parameters: as{z — 0), flxo, wi, and K. 
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FIG. 2: Evolution of /erg versus z for the constant w model with w = —0.9 and = 1. The present epoch is identified to be 
f2a:0 = 0.73 with (78(2 = 0) = 0.811. The solid line shows the numerically integrated solution, whereas the dot-dashed, dotted, 
and bold dashed lines correspond to the solutions summed up to 1-st, 3-rd, and 7-th order terms of c„ in Eq. (|52|) . 



V. VALIDITY OF ANALYTIC SOLUTIONS 



We study the validity of the analytic estimation given in the previous section. We discuss three different cases: (i) 
constant w models, (ii) tracking models, and (iii) thawing models, separately. In all the numerical simulations in this 
section, we identify the present epoch to be fi^-o = 0.73 with af,{z = 0) = 0.811. 



Constant w models 



Let us first study constant w models realized by either imperfect fluids or quintessence. Unless ^ 1 the 
dark energy perturbation is negligible relative to the matter perturbation, so that the approximate equation (j25p is 
sufficiently accurate. Since Wn =0 (n > 1), the coefficients Cn and are given, respectively, by 



(-1)" 



(54) 



and /cTg is analytically known from Eq. (|52p just by specifying the values of a^{z = 0), ^x<:>^ and w = wq. 

In Fig.[2]we plot the evolution of /erg obtained by the analytic estimation (|52p for w = —0.9 and = 1. A number 
of different lines correspond to the solutions derived by taking into account the c„ terms up to 1-st, 3-rd, and 7-th 
orders. For the growth index 7 given in Eq. (|42p we take into account the terms up to 2-nd order with respect to 
As we pick up higher-order coefficients c„, the solutions tend to approach the numerically integrated solution of 
/as. In Fig. [2] we find that the solution up to 7-th order terms of c„ can reproduce the full numerical result in good 
precision. 

In Fig. [3] wc show j a% versus z for five different values of w. In order to obtain a good convergence we pick up the 
Cn terms up to 10-th order. Figure [3] shows that our analytic estimation (|52p is sufficiently trustable to reproduce the 
numerically integrated solutions accurately. If we only pick up the terms inside 7 up to 1-st order with respect to 
rJa;, there is small difference of /erg between the analytic estimation and the numerical solutions for w > —0.6 (which 
occurs in the low redshift regime). Taking into account the 2-nd order term in Eq. (|42p . this difference gets smaller. 
Fig. [3] also displays the RSD data given in Table U which will be used to place observational constraints on w in 
Sec. Ell 
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FIG. 3: Evolution of fag versus z for the models with = 1 and (a) w — —1.2, (b) w — —1, (c) w — —0.8, (d) w = —0.6, (e) 
w = —0.4, respectively. The solid lines correspond to the numerically integrated solutions, whereas the bald dashed lines are 
derived from the analytic estimation (|52p with the 10-th order terms of c„ . We also show the current RSD data. 

0.60 I I 

0.55 - - 




0.30 - ± 

0.25 - - 

0.20 I" ■ ■ ■ I ' ' ' I ■ ' ■ I ■ ■ ' I ' ' ■ I ■ ■ ' ' 
0.2 0.4 0.6 0.8 1 12 

z 

FIG. 4: Evolution of /erg versus z for the tracking quintessence model characterized by the inverse power-law potential 
V{cj>) = M'^'^^(f>~'' . From top to bottom the solid lines correspond to the numerically integrated solutions for p = 0.1, 0.5, 1, 2, 
respectively, whereas the bold dashed lines are derived from the analytic solution H52p with the 10-th order terms. The RSD 
data are the same as those given in Fig. [S] 

B. Tracking quintessence models 

From Eq. (|33[) . we see that all coefficients w^s {n > 1) in tracking models of quintessence are expressed in terms 
of the tracker equation of state wq — il'(o)- In this case there are contributions to c„ coming from the variation of t«, 
i.e., non-zero values of /3„. Note that /3„'s depend only on W(o)- 

For the tracking quintessence with the inverse power-law potential V{(j)) — A/^+p0^p {p > 0), we compare the 
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FIG. 5: Evolution of fas versus z for the freezing quintessence with the potential V((f>) = A''[l + cos{4>/ f)]. The cases (a)-(c) 
correspond to the full numerical solutions with the model parameters (a) K — 0.81, wi = 0.17, (b) K — 1.87, wi = 4.7 x 10^^, 
and (c) K = 8.22, wi = 3.6 x 10~^, respectively, whereas the case (d) is derived by the 10-th order analytic solution (|52|l for 
K = 1.87, wi = 4.7 X 10"^ 

numerically integrated solutions of /erg with those derived by the analytic expression (|52p . In Fig. S] we show the 
evolution of /erg for p = 0.1, 0.5, 1, 2 evaluated from ([5^ as well as the numerical solutions. In Eq. ([5^ we take into 
account the c„'s up to 10-th order terms, whereas in the analytic expression of 7 in Eq. (|42|) the terms up to 2-nd 
order in fix are included. For the evaluation of fix the 1-st order solution ([55)) with flxo = 0.73 is used. From Fig. |3] 
wc find that the analytic solution ([52]) is accurate enough to reproduce the full numerical solution in high precisions. 
If we take into account the c„'s up to the 3-rd order terms, for example, there is some difference between the analytic 
and numerical results. This difference tends to disappear by including the higher-order terms of c„. While the terms 
up to 10-th order are taken into account in Fig. UJ the 7-th order solutions are sufficient for practical purpose. 

C. Thawing models 

In thawing models of quintessence the field equation of state is given by Eq. ([55]) . For larger K the deviation of 
w from —1 occurs at later times with a sharper transition. From Eq. (j37p we find that the higher-order terms in 
are not negligible for K larger than the order of 1. In fact we have numerically confirmed that, for K > 0(1), the 
expansion (|37p docs not accommodate the rapid transition of w at late times unless the higher-order terms arc fully 
taken into account. This property also holds for the evolution of fa^. Only when K is smaller than the order of 1, 
the analytic estimation (j52p can reproduce the full numerical solutions in good precision. 

In Fig. [5] we plot the numerical evolution of /erg for the potential V{(f>) — A^[l + cos(0//)] with three different 
values of K. Since w is close to —1 in all these cases and the deviation from w = —1 occurs only at late times, the 
evolution of /erg is not very different from each other for K < 10. From this analysis, it is clear that, only when very 
accurate data of /tjg are available in the future, it will be possible to distinguish between the models with different 
values of K. In Fig.[5]we also show the analytic solution derived for K = 1.87 and wi = 4.7 x 10~^ as the bold dashed 
line (d). We take into account the c„'s up to 10-th order to evaluate /erg in Eq. ([5^ . Compared to the full numerical 
solution labelled as (b), there is a small difference in the high-redshift regime. We confirm that this deviation tends 
to be smaller by involving the c„'s higher than 10-th order. For K < 1 the analytic estimation is more accurate even 
without including such higher-order terms. 
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VI. CONSTRAINTS FROM THE CURRENT RSD DATA 



111 this section, we place observational bounds on two models of dark energy discussed in Sees. IV Al and IV Bl by 
using the current RSD data presented in Table HI For the today's value of erg wc consider the prior obtained from 
observations of CMB, BAO, and Hubble constant measurement (Hq), i.e., 

cr8(z = 0) = 0.816 ±0.024. (55) 

Here and in what follows all the error bars correspond to the 68.3% confidence level (CL). Recall that we derived the 
analytic formula (|52p under the approximation that the dark energy perturbation is neglected relative to the matter 
perturbation. For the validity of this approximation we put the prior w < —0.1. 



A. Constant w models 



For the models of constant w the today's matter density parameter constrained from SN la, CMB, BAO, and Hq 
observations is 0] 

Qmo = 0.272±° °i3 . (56) 
Under the priors (|55|) and ([56)) we estimate the best-fit to the set of parameters P = {w, a%(z — 0), il,no) by evaluating 

2 /o 

the likelihood distribution function, L oc e^^ ' , with 

2 /./cr8,ob(Zj) - ./0-8,th(^^i)\^ ,r-7\ 

^ ^gl ) • 

Here f(Js.oh{zi) {i = 1, ■ ■ ■ ,9) are the 9 data displayed in Table H] with the error bars (Ti, whereas f(Ts,th{zi) are the 
theoretical values derived from the analytic solution (|52p. For the evaluation of fcrg thizi) we pick up the c„'s up to 
10-th order. For the growth index 7 the terms up to 2-nd order with respect to fix are included in Eq. (j42|). For 
fix we use the 1-st order solution (|53l). In tracking quintessence models analyzed in Sec. IVIBI we also take the same 
orders of expansions for /ctsi 7i a-nd fix- 

We find that the best-fit model parameters are 

w = -0.604, (78(2 = 0) = 0.840, r2„o = 0.285, (58) 

with reduced Xr = 0.947 {Xr = Xmin/'^' where ly stands for the degrees of freedom). At 68.3% CL, our analysis 
restricts the equation-of-state parameter to the interval 

- 1.245 < itx -0.347, (59) 

whereas as{z = 0) and flmo are unconstrained by current data even assuming the priors ()55P and ()56|) . Although 
the current bounds on w are weaker than those arising from background tests (see, e.g., Refs. 0,0]), we expect that 
upcoming RSD data from ongoing and planned galaxy redshift surveys can improve this situation in the near future. 



B. Tracking quintessence models 



For the tracking quintessence models in which the equation of state is given by Eq. ([5^ we also carry out the 
likelihood analysis by using the analytic solution ([52)) of fas as well as the expressions for 7 and ilx given in Eqs. ([42)) 
and ()53p respectively. While the equation of state ([5^ is derived for quintessence, we do not impose the prior 
""^(o) > ~1 for generality. For this kind of models, a joint analysis involving current SN la, CMB, and BAO gives the 



following bound on the matter density parameter |6l| : 



0.273 < a^o < 0.293 , (60) 
which is used in our analysis as a prior for r2,„o. The best- fit model parameters are found to be 

u;(0) = -0.461 , (78(2 = 0) = 0.840, r2™o = 0.293 , (61) 



13 



with xl = 0.923. At 68.3% CL, we found 

- 1.288 < W(o) < -0.214, (62) 

whereas the parameters cFsi^ = 0) and flmo are again unconstrained in the regions of (|55p and (|60|). As expected, 
the bounds on ii;(Q) are weaker than those obtained in constant w models [Eq. (|59p ]. We note that in tracking models 
the equation of state decreases at late times, which is accompanied by the decrease of fag. Compared to constant w 
models, this allows the possibility to fit the data better even for larger values of w during the matter era. 

VII. CONCLUSIONS 

In this paper we have provided an analytic formula of fas for dynamical dark energy models in the framework of 
GR. This was derived by using the approximate matter perturbation equation (j24p . which is trustable as long as the 
contribution of the dark energy perturbation to the gravitational potential is negligible relative to that of the matter 
perturbation. Our formula of /erg can be applied to many dark energy models including imperfect fluids, quintessence, 
and k-essence in which the sound speed squared Cg is not very close to 0. 

Our derivation of /erg is based on the expansion of w with respect to the dark energy density parameter fl^, i.e., 
w ^ wq + J2n=i ''^ni^x)"- The growth rate / = S'^/Sm of matter perturbations is parametrized by the growth index 
7, as / = (1 — ilx)"'- We expanded 7 in terms of fix up to 2-nd order terms. Since 7 is dominated by the term 
7o = 3(1 — wo)/(5 — 6wo), it is a good approximation to treat 7 as a constant for the derivation of the integrated 
solution of dm- The c„'s in Eq. (|52p are given by Eq. (|47p . where a„ and /3„ appear as the coefficients of the expansion 
of the terms (1 — and 1/w respectively. For the density parameter fl^, the 1-st order solution ((53)) is usually 

sufficient to obtain accurate analytic solutions of fag. 

In Sec. |V]we have studied the validity of the analytic formula (|52|) in concrete models of dark energy. For constant 
w models in which c„ and fix are given by Eq. (|54p . the analytic solution up to 7-th order terms of c„ reproduces 
the numerically integrated solutions with good precision. This property also holds for tracking quintessence models 
where the evolution of w is given by Eq. (|33l). In thawing quintessence and k-essence models, where w is given by 
Eq. (|35l). the formula (|52p can be trustable K < 1, but for K larger than the order of 1, we need to fully take into 
account the higher-order terms of c„ to have good convergence of fas- 

In Sec. lVII we have discussed observational constraints on two different dark energy models by using the current RSD 
data. In both constant w and tracking quintessence models the analytic solution (|52p includes the three parameters 
as{z = 0), fimOi and w (or u)(o)). Under the priors on (7s{z = 0) and flmo constrained by SN la, CMB, BAO, and 
Hq measurements, we derived the bounds —1.245 < w < —0.347 (68%CL) for constant w models and —1.288 < 
10(0) < —0.214 (68%CL) for tracking quintessence models. Although the upper bounds on the dark energy equation 
of state are still weak with current data, we expect to obtain more precise data from ongoing surveys or near-future 
projects such as Subaru/FMOS, HETDEX, and J-PAS. We hope that our analytic formula of fas will be useful to 
place tighter bounds on dynamical dark energy models in the future. 
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